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Motivated by some recent criticisms to our alternative Langevin equation for driven lattice gases 
(DLG) under an infinitely large driving field, we revisit the derivation of such an equation, and 
test its validity. As a result, an additional term, coming from a careful consideration of entropic 
contributions, is added to the equation. This term heals all the recently reported generic infrared 
singularities. The emerging equation is then identical to that describing randomly driven diffusive 
systems. This fact confirms our claim that the infinite driving limit is singular, and that the main 
relevant ingredient determining the critical behavior of the DLG in this limit is the anisotropy and 
not the presence of a current. Different aspects of our picture are discussed, and it is concluded that 
it constitutes a very plausible scenario to rationalize the critical behavior of the DLG and variants 
of it. 
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The driven lattice gas (DLG) |l]-|j| is a simple non- 
trivial extension of the kinetic Ising model, and consti- 
tutes certainly a main archetype of out-of-equilibrium 
system. Fully understanding the critical properties of the 
DLG would be a fundamental milestone on the way to 
rationalizing the fast developing field of nonequilibrium 
phase transitions. The DLG is defined as a half filled, d- 
dimensional kinetic Ising model with conserved dynam- 
ics, in which transitions in the direction (against the di- 
rection) of an external field, E, are favored (unfavored) 
. The external field induces two main nonequilib- 
rium effects: the presence of a net current of particles in 
its direction, and anisotropic system configurations. At 
high temperatures the system is in a disordered phase, 
while below a certain critical point it orders by segre- 
gating into high and low density aligned-with-the-field 
stripes. 

In order to analyze the DLG critical nature, and de- 
termine its degree of universality, a Langevin equation 
intended to capture the relevant physics at criticality was 
proposed and renormalized more than a decade ago [Q. 
This elegant theory, the driven diffusive system (DDS) 
seems to capture the main symmetries and conservation 
laws of the discrete DLG (including a current term as 
the most relevant nonlinearity) , and is therefore a suit- 
able and very reasonable candidate to be the canonical 
continuous model, representative of the DLG universality 
class. 

Unfortunately, the most emblematic prediction com- 
ing from the analysis of the DDS equation, namely, the 
mean field behavior of the order parameter critical ex- 
ponent (/?= 1/2 ||), has not been compellingly verified 
in any Monte Carlo simulation of the DLG in spite of 
the huge computational effort devoted to test it. In par- 
ticular, systematic deviations from scaling are observed 
both in d = 2 |||| and in d = 3 M if data collapse is 
attempted using = 1/2 0. On the other hand, dif- 



ferent Monte Carlo numerical simulations (performed in 
different geometries and using different finite size scal- 
ing ansatzs) lead systematically to a value of around 
0.3 with error bars apparently excluding — 1/2 (we 
refer the interested reader to || for a review of simula- 
tion analysis). This is a main indication that the DDS 
Langevin equation does not describe properly the DLG 
at criticality. 

Moreover, there are some other hints suggesting 
strongly that the discrepancies between the predictions 
of the standard theory and Monte Carlo results are more 
fundamental than a simple numerical difference in 0. 
In particular, the intuition developed from Monte Carlo 
simulations of the DLG and variants of it |||| (performed 
under large external driving fields) suggests that, con- 
trarily to what happens for the DDS equation, it is the 
anisotropy and not the presence of a current the most 
relevant ingredient for criticality. For instance, in a mod- 
ification of the DLG in which anisotropy is included by 
means others than a current ||, the scaling behavior at 
criticality remains unaltered upon the switching on of an 
infinite driving (see the appendix and ||[|). Other com- 
pelling evidences supporting this hypothesis can be found 

in IJnJ. 

In order to shed some light on this puzzling situation 
and reconcile theory with numerics, different possible sce- 
narios have been explored; but so far no satisfactory clar- 
ification has been reached. Within this context, we have 
recently revisited the time-honored DDS equation and 
questioned its general validity jll],[l2). In particular, we 
have tackled the task of constructing a coarse-grained 
procedure in a more detailed way such that, starting from 
a Master equation representing the DLG, would give as 
output a continuous Langevin equation. This approach 
permits us to keep track of microscopic details that could 
eventually be overlooked when writing down a Langevin 
equation respecting naively the microscopic symmetries 
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and conservation constraints m . This approach has given 
rise to a rather unexpected and quite interesting out- 
put: The limit of infinitely large driving (i.e. the limit 
in which attempted jumps in the direction of the field 
are performed with probability one and jumps against E 
are strictly forbidden) is singular Let us stress that 
in order to enhance nonequilibrium effects most of the 
available computer studies are performed in this limit. 
The main results derived so far using our approach are: 

• For vanishing values of the driving field it leads to 
the standard equilibrium model B p3| , capturing 
the relevant physics of the kinetic Ising model with 
conserved dynamics. 

• For nonvanishing, but finite driving fields we re- 
produce the standard DDS Langevin equation 

• In the limit of infinitely large driving, where the de- 
pendence of jumps in the direction of the field on 
energetics is replaced by a zero-one (all-or-nothing) 
condition, a different Langevin equation emerges. 
This new equation has the main property of not 
including any relevant term coupling E to the den- 
sity field <p Ml ; an d the presence of anisotropy is 
its main relevant ingredient. 

The new Langevin equation for the DLG under in- 
finitely large driving proposed in M ; P^ 1 an d renormal- 
ized in has some important virtues to be discussed 
afterwards, but seems also to exhibit some pathologies, as 
recently pointed out by Caracciolo et al. |l5| and also by 
Schmittmann et al. |l6| |. In what follows we show that 
such anomalies can be healed in a rather natural way, 
and do not disprove at all the general validity of our new 
approach (as could be inferred from (l(|). 

Let us now present the Langevin equation derived in 
|l|>@ for the infinite dr iving limit, report on its defi- 
ciencies |lf|[l6|], and discuss the way to heal them. The 
equation reads |TI|Jl^]: 



the density at each column. These (spurious) conserva- 
tion laws, absent in the DLG, are at the origin of the 
infrared singularities appearing in Eq. (|l]) ]l&|L{l. 

In order to investigate the causes of this deficiency in 
our Langevin equation and eventually overcome the prob- 
lem of the extra conservation laws and associated infrared 
divergences, we have re-analyzed our derivation of Eq. 
(0) in M i B ■ O ne can easily see that the transition rates 
in the microscopic master equation in pT| , |l2"| ] were written 
as depending on the variations of two adding contribu- 
tions: the free energy functional (the usual Ginzsburg- 
Landau free energy) and the external driving-field con- 
tribution. The transition rates, written in that way, sat- 
urate to zero or one in the field direction, in the limit 
of infinite driving. This saturation erases any further de- 
pendence on the free energy density (which includes both 
entropic and energetic contributions). On the contrary, 
in the DLG it is only the dependence on the Ising ener- 
getics that becomes negligible in the limit of large driving 
fields. In a coarse grained description we should therefore 
separate energetic from entropic terms. With this guid- 
ing idea, we have reconsidered our derivation of Eq. (|l|) 
and rewritten the transition rates in [jll 12 as the prod- 
uct of two contributions: one controlling the energetics 
and the other one the entropic par t WA- By performing 
a calculation analogous to that in [[Ll| , |12| , but including 
the transition rates written in this modified way, it is 
a matter of algebra to see that a new term (missing in 
@PHl!) emerges: pV||0(x,t) @. 

It is straightforward to verify that apart from prop- 
erly keeping track of entropic contributions, this extra 
(mass) term heals all of the aforementioned problems in 
Eq. (Jy): no spurious conservation laws are involved and 
generic infrared singularities disappear. 

Let us now discuss how this new additional term affects 
the results presented in [jlJJ . Performing a naive scaling 
analysis, one sees that xn ~ x± 2 |l9| , and upon elimi- 
nation of naively irrelevant terms and absorbing eo into 
the time scale, one obtains our final result: the critical 
Langevin theory under infinitely large driving 
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where Vn (Vx) is the gradient operator in the direction 
parallel (perpendicular) to the driving field, and £ is a 
conserved Gaussian white noise [pr| , pT[ . This equation 
is analogous to a model B in the direction(s) perpendic- 
ular to the field, coupled to a simple random diffusion 
mechanism in the parallel direction. The origin of all 
the difficulties pointed out in [ p^|JT^ ] can be traced back 
to the following property: Defining the total density for 
each value of ry, T(ru,t) = J d d ~ 1 r±cf)(r\\,r±), it is not 
difficult to see (after averaging over the noise) that T(ry) 
is a conserved quantity for all values of ry | [l5fl . Observe 
also that T(ry) is nothing but the zero Fourier mode of 
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that we call the anisotropic diffusive system (ADS). This 
turns out to be a well known Langevin equation: the 
continuous representation of the randomly driven DLG 
pojpl], i.e. a DLG in which the external field changes sign 
randomly in an unbiased fashion. The main difference 
between this theory and the DDS is that the ADS does 
not include an overall current. The current term i?V||0 2 
appearing in the DDS (and constituting its most relevant 
nonlinearity) is absent here. In the random DLG, such a 
term cannot appear for symmetry reasons, while in the 
infinite driving case discussed in this paper it is the sat- 
uration of the transition rates in the field direction that 
prevents such a current term from appearing. 
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The cubic operator and the Laplacian term in the par- 
allel direction in Eq. (0) are both marginal at the critical 
dimension d = 3. The results up to first order in an ep- 
silon expansion of Eq. (Q) around d = 3 are po| , pl| : 
vx = 1/2 + E/T2, and [3 = 1/2 - e/6 Observe 
that in cZ = 2 one obtains /? = 1/3 (slightly modified 
by two-loop corrections pof ) in remarkable good agree- 
ment with Monte Carlo results. For instance (see table 
and p3[): the best available Monte Carlo result for the 
random DLG is (3 w 0.33 f20|] ; for the infinitely driven 
DLG [3 » 0.30 ± 0.05 @; and /? » 0.34 for the closely 
related model studied in ||, called ALGA, and argued 
to belong to the same universality class (see appendix). 

Some further comments on the validity of our approach 
follow: (i) Our complete theory (including constant and 
irrelevant terms) does have a net current pl| , [T2] | , though 
it does not enter the final Langevin equation, (ii) An 
infinitely large field is, in practice, any for which transi- 
tions against the field never occur. Given that all com- 
monly used transition rates depend on E through ex- 
ponential functions, field values much larger that unity 
can be considered infinite for all practical purposes in 
Monte Carlo experiments. For smaller fields, we expect 
crossover effects from the infinite field regime (ruled by 
the ADS) to the finite-driving standard DDS behavior 
to occur. These crossovers could obscure the numeri- 
cal observation of the DDS mean-field exponent beta for 
large but finite driving fields, (iii) The introduction of 
the new term in the parallel direction heals all the pos- 
sible problems in relation to infrared singularities, extra 
conservation laws, and anomalies in the structure func- 
tion ||l6| . In particular, the structure function presents a 
discontinuity singularity as happens in the DLG g] . (iy) 
Given the absence of any relevant current term in Eq. (g), 
the critical theory has "up-down" symmetry (0 — > —<j)). 
This symmetry, in principle, is absent in the microscopic 
model as the presence of nonvanishing three-point cor- 
relation functions seems to indicate B]. However, it is 
not clear whether such correlations are relevant at crit- 
icality or not. As an indication that in fact they could 
well be irrelevant, we discuss here the problem of trian- 
gular anisotropies ^2|: In both the DLG and the DDS, 
droplets of the minority phase (if any) develop triangular 
shapes, closely related to the existence of nonvanishing 
three-point correlations. However, triangles orientate in 
opposite directions in the microscopic DLG and in the 
continuous DDS. This difference seems to be not univer- 
sal as shown by recent Monte Carlo studies j22], i.e. it 
depends on microscopic details and can be modified by 
changing them both in the DLG and in the DDS. This 
fact supports the idea that nonvanishing three-point cor- 
relation functions might not be a relevant ingredient for a 
description of the DLG at criticality. More significatively, 
simulations show that for large enough driving fields the 
triangular anisotropy is suppressed (see 22 1), providing 



an indication that the up-down symmetry is restored in 
the infinite driving limit. This constitutes, we believe, 
another strong backing of our picture. 



In summary, we have discussed the plausibility of the 
alternative field theoretical approaches to driven lattice 
gases under the effect of an infinitely large external driv- 
ing field. Some deficiencies recently pointed out are over- 
come by introducing an extra Laplacian term in the direc- 
tion of the field in the Langevin equation first proposed 
in pd|Jl^ ]. This new term, coming from a proper con- 
sideration of entropic contributions, had been overlooked 
in previous papers. Our approach leads to the follow- 
ing global picture: (i) For E = 0, model B reproduces 
the equilibrium critical properties of isotropic diffusive 
systems, (ii) For finite driving field, the standard DDS 
Langevin equation, including a current term, should de- 
scribe properly the long wavelength properties around 
the critical point, (iii) The limit of infinitely fast driving 
is singular: the current is irrelevant, and the anisotropy 
becomes the main relevant property. The leading critical 
properties in this case are expected to be described by 
the ADS, Eq. (||). The reason for this being that in the 
presence of infinite driving the transition rates saturate 
to 1 (0) for allowed (forbidden) transitions in the driv- 
ing direction, and no further track of coupling between E 
and the density field survives in the resulting Langevin 
equation. We want to stress at this point that this prop- 
erty was not obvious a priori, but emerges as a natural 
output from our model building strategy. 

The proposed Langevin equation for large external 
field provides a quite plausible scenario shedding some 
light on a difficult problem. In particular, it justifies the 
observed lack of differences (for large fields) in simula- 
tions in systems with and without a current, and pro- 
vides a likely justification of why the standard prediction 
(3 = 1/2 is not confirmed in Monte Carlo simulations, 
and instead a value [3 » 0.33 is observed. 

In order to test numerically the picture presented in 
this paper, it would be highly desirable to perform ex- 
tensive simulations for finite driving field (E w 1), and 
study whether differences with respect to the available 
Monte Carlo results for large fields appear. It would also 
be interesting to improve the finite size scaling analysis 
following the strategy used in |24|,[l5| . 

Appendix. 

As an evidence aimed to transmitting the intuition 
that, for infinitely large driving, the current is not rel- 
evant at criticality, let us briefly discuss in this ap- 
pendix a rather compelling Monte Carlo observation. It 
corresponds to a variation of the DLG, named ALGA 
(anisotropic lattice gas automaton); see || for a detailed 
definition. This model is placed by definition at the limit 
of infinite driving: jumps in the anisotropy direction are 
performed randomly without attending to energetic con- 
siderations. Simulations are performed both in the pres- 
ence of an overall current (case p ^ 1/2 in ||) and in 
the absence of it (p = 1/2); the curves for the order 
parameter versus the distance to the critical point are 
indistinguishable in the cases with and without a current 
(figure 3 in (9) is particularly illuminating). It could be 
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argued that the details of this modified model M render 
it not completely equivalent to the original DLG. How- 
ever, we do not think these microscopic differences have 
any relevance at a coarse grained level. In fact, we expect 
this model to be represented by Eq. (^) : in one direction 
particles tend to stay together, and it is natural to as- 
sume that their coarse grained behavior is controlled by 
a model B in this direction. In the other direction, jumps 
occur regardless of energetics and, therefore, the dynam- 
ics becomes purely diffusive. With these two ingredients 
we recover the ADS, Eq. (^]), as the Langevin equation 
for the ALGA. As a further evidence supporting this hy- 
pothesis let us mention that the measured (3 exponent 
in the ALGA is (3 w 0.34 (again very close to the value 
(3 1/3) in both cases: with and without current. 
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MODEL 




No net current 


Net current 
Finite driving 


Net current 
Infinite driv. 


DLG 




Model B 
/3 = 0.125 


DDS 


ADS 
f3 w 0.3 


Random DLG 




ADS 
P w 0.33 


DDS 


ADS 


ALGA 




ADS 
S3 w 0.34 


Undefined 


ADS 
/? » 0.34 



Universality classes for different models as a function of the 
driving field value. The reported values of f3 correspond to 
the best-to-date numerical results in d = 2. The theoretically 
predicted value is j3 = 1/2 for the DDS and j3 ft! 0.33 for the 
ADS. 
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